Received: 7 November 2016 


DOI: 10.1111/tbed. 12630 

ORIGINAL ARTICLE 


Wiley 


Transboundary and Emerging Diseases 

T7 


Coronavirus infections in horses in Saudi Arabia and Oman 

M. G. Hemida 1,2+ | D. K. W. Chu 3+ & | R. A. P. M. Perera 3+ | R. L. W. Ko 3 | 

R. T. Y. So 3 | B. C. Y. Ng 3 | S. M. S. Chan 3 | S. Chu 3 | A. A. Alnaeem 4 | 

M. A. Alhammadi 1 | R. J. Webby 5 | L. L. M. Poon 3 | U. B. R. Balasuriya 6,7 | M. Peiris 3 


department of Microbiology and 
Parasitology, College of Veterinary Medicine, 
King Faisal University, Al-Hasa, Saudi Arabia 

department of Virology Faculty of 
Veterinary Medicine, Kaferelsheik University, 
Kaferelsheik, Egypt 

3 School of Public Health, The University of 
Hong Kong, Hong Kong, China 

department of Clinical Studies, College of 
Veterinary Medicine, King Faisal University, 
Al-Hasa, Saudi Arabia 

division of Virology, Department of 
Infectious Diseases, St. Jude Children’s 
Research Hospital, Memphis, TN, USA 

6 Maxwell H. Gluck Equine Research Center, 
Department of Veterinary Science, College 
of Agriculture, Food and Environment, 
University of Kentucky, Lexington, KY, USA 

department of Microbiology, Immunology 
and Molecular Genetics, College of Medicine, 
University of Kentucky, Lexington, KY, USA 

Correspondence 

M. Peiris, School of Public health, The 
University of Hong Kong, Pokfulam, Hong 
Kong Special Administrative Region, China. 
Email: malik@hku.hk 

Funding information 

King Abdulaziz City for Science and Technol¬ 
ogy, Grant/Award Number: ARP-34-117; US 
National Institute of Allergy and Infectious 
Diseases, Grant/Award Number: 
HHSN272201500006C; National Institutes 
of Health, Department of Health and Human 
Services; Health and Medical Research Fund; 
Food and Health Bureau, Government of the 
Hong Kong Special Administrative Region; 
King Abdul-Aziz city of Science and Technol¬ 
ogy (KACST), Grant/Award Number: ARP- 
34-117 


Summary 

Equine coronaviruses (ECoV) are the only coronavirus known to infect horses. So 
far, data on ECoV infection in horses remain limited to the USA, France and Japan 
and its geographic distribution is not well understood. We carried out RT-PCR on 
306 nasal and 315 rectal swabs and tested 243 sera for antibodies to detect coron¬ 
avirus infections in apparently healthy horses in Saudi Arabia and Oman. We docu¬ 
ment evidence of infection with ECoV and HKU23 coronavirus by RT-PCR. There 
was no conclusive evidence of Middle East respiratory syndrome coronavirus infec¬ 
tion in horses. Serological data suggest that lineage A betacoronavirus infections are 
commonly infecting horses in Saudi Arabia and Oman but antibody cross-reactivities 
between these viruses do not permit us to use serological data alone to identify 
which coronaviruses are causing these infections. 
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INTRODUCTION 


Coronaviruses infect a wide range of mammals and birds, exhibit a 
marked tropism for epithelial cells of the respiratory and enteric 


iThese authors contributed equally to this manuscript and are co-first authors. 


tracts and cause a remarkably diverse spectrum of diseases in 
humans and livestock. They typically have a restricted host range, 
infecting only their natural host and closely related animal species 
but do have the capacity to cross the species barrier to infect new 
hosts. Equine coronavirus (ECoV) is the only coronavirus known to 
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infect or cause disease in horses (Balasuriya, 2013). It belongs to the 
genus Betacoronavirus lineage A, as does human coronavirus OC43, 
bovine coronavirus and porcine hemagglutinating encephalomyelitis 
virus. Equine coronavirus was first isolated from faeces of a diar- 
rhoeic foal in 1999 (ECoV-NC99) in North Carolina, USA (Guy, Bres- 
lin, Breuhaus, Vivrette, & Smith, 2000), and was initially believed to 
only affect foals. Since 2010, there have been several reports of 
ECoV-associated respiratory and enteric infections in adult horses in 
Japan, Europe and the United States, but its global distribution is still 
poorly defined (Kooijman, Mapes, & Pusterla, 2016; Miszczak et al., 
2014; Oue, Morita, Kondo, & Nemoto, 2013; Pusterla, Holzen- 
kaempfer, Mapes, & Kass, 2015). 

Middle East respiratory syndrome coronavirus (MERS-CoV) is a 
betacoronavirus within lineage C first identified in humans in 2012 
and continues to pose a threat to global health (WHO, 2016 May 
27). The evidence points to dromedary camels as a natural host 
and a major source for human MERS-CoV infections (Alexandersen, 
Kobinger, Soule, & Wernery, 2014; Chu et al., 2015; Gutierrez, 
Tejedor-Junco, Gonzalez, Lattwein, & Renneker, 2015; Memish 
et al., 2014; Perera et al., 2013; Reusken et al., 2014). There is no 
convincing evidence of MERS-CoV infections in other domestic 
livestock species so far (Adney et al., 2016; Alexandersen et al., 
2014; Hemida et al., 2013; Meyer et al., 2015; Perera et al., 2013). 
Sequence similarity comparisons of dipeptidyl peptidase-4 (DPP4 
[CD26]), the functional receptor for MERS-CoV, have revealed that 
equine DPP4 is phylogentically closely related to human DPP4 and 
the binding affinity of MERS-CoV spike SI domain for equine 
DPP4 is similar to that of human and camel DPP4, raising the pos¬ 
sibility that MERS-CoV infections may occur among horses (Barlan 
et al., 2014). Serological studies of horses originating from Spain 
and the United Arab Emirates gave negative results for MERS-CoV 
(Alexandersen et al., 2014; Meyer et al., 2015), but it is still rele¬ 
vant to investigate for evidence of MERS-CoV in horses in areas 
endemic for MERS-CoV. 

We have tested for coronavirus infections in apparently 
healthy horses in Saudi Arabia and Oman by testing nasal and 
rectal swabs by RT-PCR assays for MERS-coronavirus and using a 
pan-coronavirus RT-PCR with potential to detect all coronaviruses, 
hitherto known and unknown. Following the detection of ECoV 
and HKU23 coronavirus (HKU23) by the pan-coronavirus RT-PCR, 
we tested the swab specimens with specific RT-PCR assays for 
ECoV and HKU23 and tested the sera with serological assays to 
detect ECoV, bovine CoV (BCoV) (closely related to HKU23), as 
well as MERS-CoV. 


2 


MATERIALS AND METHODS 


2.1 | Specimen collection 


Serum samples (n = 243), nasal swabs (n = 306) and rectal swabs 
(n = 315) were collected from apparently healthy adult (>1 years of 
age) horses from East and Central provinces of the Kingdom of 
Saudi Arabia and Oman. Many horses had all three samples 


collected. These specimens were collected during the period Febru¬ 
ary 2014 to January 2015. The numbers of specimens collected from 
different geographic regions of Saudi Arabia and Oman are indicated 
in Table 1. The animal management practices and exposures to other 
domestic animal species such as dromedaries are indicated. Farming 
activities include activities related to date palm plantations as well as 
other farming activities. Racing activities refer to horse racing within 
Saudi Arabia rather than participation in international races. 


2.2 | Cells and viruses 

Vero cells (ATCC CCL-81) were used to culture and propagate 
MERS-CoV and the HRT-18G cell line (ATCC, CRL-11663, American 
Type Culture Collection, Manassas, VA, USA) for culture of ECoV 
and BCoV. The culture medium for both the cell lines was DMEM 
(Cat:10569; Life Technologies, Carlsbad, CA, USA) supplemented 
with 10% foetal calf serum (Cat:16140; Life Technologies). 

The prototype viruses used were MERS-CoV (EMC/2012 
obtained from Dr R Fouchier, Erasmus MC, Rotterdam), the ECoV- 
NC99 strain (obtained from Dr. Udeni B.R. Balasuriya, Maxwell H. 
Gluck Equine Research Center, University of Kentucky, USA) (Guy 
et al., 2000; Zhang et al., 2007) and the bovine coronavirus (ATCC 
BRCV-OK-0514-2). 

2.3 | RT-PCR testing 

Viral RNA was extracted from nasal and rectal swabs using the 
QIAamp viral RNA minikit (Qiagen, Hilden, Germany) following the 
manufacturer’s instructions. RNA extracts were tested for evidence 
of conserved coronavirus nucleic acid genetic sequences using previ¬ 
ously reported RT-PCR assays (Chu et al., 2014), RTqPCR assay for 
MERS-CoV upE gene (Corman et al., 2012), RTqPCR assay for ECoV 
(Miszczak et al., 2014), and a RTqPCR assay for HKU23 reported 
below. RNA was reverse transcribed in a 20 pi reaction mixture con¬ 
taining lx first-strand buffer, 5 mM DTT, 0.5 mM deoxynucleotide 
triphosphates (dNTP), 2.5 ng/pl random hexamers and 200 units of 
Superscript III (Life Technologies). 

The pan-coronavirus nested PCR targeted the RNA-dependent 
RNA polymerase (RdRp) gene of coronaviruses developed by us as 
previously described (Chu et al., 2014). Using the cDNA synthesized 
as described above, first round PCR was carried out using forward 
primer 5'-GGKTGGGAYTAYCCKAARTG-3' (position 15,287 of ECoV 
strain NC99; GenBank accession number EF446615) and reverse pri¬ 
mer 5'-TGYTGTSWRCARAAYTCRTG-3' (position 15,869 of ECoV 
strain NC99) and Phusion High Fidelity PCR Master Mix Kit (Thermo 
Scientific). A 25 pi reaction mixture contained 12.5 pi of 2x phusion 
master mix, 1.5 pi (10 pmol) forward primer, 1.5 pi (10 pmol) reverse 
primer, 4.5 pi ddH 2 0 and 5 pi of cDNA. The PCR cycler conditions 
for the amplification were 98°C for 2 min, 40 cycles of amplification 
(98°C for 15 s, 48°C for 15 s, 72°C for 30 s), then 72°C for 2 min. 
An aliquot of 0.5 pi of the first PCR product was then further ampli¬ 
fied in a second round PCR prepared as first around PCR but using a 
new set of primers (forward primer 5'-GGTTGGGACTATCC 
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TABLE 1 Geographic location of specimens collected and animal management practices 


Region 

City 

Provider 

code 

Animal 

management 

practices 

Contact with 

other animals 

Animals in contact with 

Nasal 

swabs 

Rectal 

swabs 

Sera 

Virus RNA 
detected, if any 

Saudi Arabia 
(Central) 

Riyadh 

1 

Breeding 

Occasionally 

Camels, sheep, goat, chickens 

14 

14 

14 

ECoV (n = 5) 



2 

Farming, 

Breeding 

Occasionally 

Camels, sheep, goat 

0 

0 

33 


Saudi Arabia 
(East) 

Dammam 

3 

Farming, 

Breeding 

Frequent 

Camels, sheep, goat, chickens 

103 

103 

19 

HKU23 



4 

Breeding 

Frequent 

Camels, sheep, goat 

6 

6 

6 


Saudi Arabia 
(East) 

Al-Hasa 

5 

Farming, 

Breeding 

Occasionally 

Camels, sheep, goat 

39 

39 

39 




6 

Farming 

Frequently 

Camels, sheep, goat, chickens 

27 

27 

0 




7 

Farming 

Frequently 

Camels, sheep, goat, chickens, dogs 

26 

26 

15 




8 

Farming, 

Breeding 

Frequent 

Camels, sheep, goat, chickens 

22 

22 

12 


Saudi Arabia 
(Central) 

Al-Qassem 

9 

Racing 

Rare 

Camels, sheep, goat 

20 

19 

19 


Saudi Arabia 
(East) 

Qateef 

10 

Farming 

Frequent 

Camels, sheep, goat, chickens 

41 

41 

24 

HKU23 



11 

Farming 

Frequent 

Camels, sheep, goat, chickens 

0 

0 

13 


Saudi Arabia 
(East) 

Aljobail 

12 

Farming 

Frequent 

Camels, sheep, goat, chickens 

0 

0 

5 


Oman 


13 

Racing, 

breeding 

Frequent 

Camels, dogs 

8 

18 

44 


Total 






306 

315 

243 



TAAGTGTGA-3' (position 15,287 of ECoV strain NC99), reverse pri¬ 
mer 5'-CCATCATCAGATAGAATCATCAT-3' (position 15,721 of 
ECoV strain NC99) which amplify a final PCR product of 440 bp. 
The final PCR amplicons (440 bp) were gel purified using QIAquick 
gel purification kit (QIAGEN) and analysed by DNA sequencing using 
the primers used for the second round PCR. To amplify an extended 
RdRp gene sequence of both equine CoV and HKU23-like viruses, 
PCR was prepared as above with forward primer 5'-TACGCAG 
GATGGTAATGCTG-3' (position 14,767 of ECoV strain NC99) and 
reverse primer 5'-TTGTGCGCATTCATTCGCAAG-3' (position 15,404 
of ECoV strain NC99) and 5 pi of cDNA as template for the reac¬ 
tion. PCR products with expected size of 658 bp were purified for 
DNA sequencing using the forward and reverse primers. 

We developed a two-step real-time quantitative RT-PCR (RT- 
qPCR) assay for the detecting HKU23 targeting N gene of the virus. 
Virus RNA was reverse transcribed as detailed above. Real-time PCR 
was performed using Power SYBR Green Master Mix (Life Technolo¬ 
gies) in a 20 pi reaction containing 5 pmol forward primer 5'- 
CTGCCACGATGGTATTTTTAC-3' (position 29,758 of HKU23 strain 
362F) and 5 pmol reverse primer 5'-GGTATTGACAT 
CAGCCTGGT-3' (position 29,852 of HKU23 strain 362F) with 40 
cycles of amplification (95°C for 15 s, 62°C for 40 s) in a 20 pi vol¬ 
ume followed by melting curve analysis. Positive result was inter¬ 
preted by a detectable Ct value >40 considered together with a 
amplicon melting temperature within the range of 81.64°C to 
82.64°C. 


The details of all RT-PCR and real-time PCR assays used in this 
study are summarized in Table SI. 


2.4 


PCR amplicon sequencing and analysis 


The primers used for sequencing were those used to generate the 
PCR products. DNA amplicons were purified and sequenced in both 
forward and reverse directions with the PCR primers using BigDye 
Terminator v3.1 Cycle Sequencing Kit (Applied Biosystems, Foster 
City, CA) and a 96-capillary 3730x1 DNA Analyzer (Applied Biosys¬ 
tems) according to manufacturer’s instructions. Sequences were 
assembled using BioEdit 7.0 and MUSCLE for multiple sequence 
alignment. Phylogeny analysis was performed by maximum likelihood 
method using PhyML with GTR nucleotide substitution model and 
SH-like approximate likelihood ratio test for estimating branch 
support. 


2.5 | Serology 

All sera were tested using previously validated MERS-CoV spike pseu¬ 
doparticle neutralization test (ppNT) (Hemida et al., 2014; Perera 
et al., 2013), and the positive sera were confirmed using microneutral¬ 
ization (MN) test and plaque reduction neutralization tests (PRNT) 
using live MERS-CoV strain EMC in a Biosafety Level 3 containment 
laboratory (Park et al., 2015). Selected sera were tested for antibodies 
to ECoV and BCoVs (Perera et al., 2013). The North American ECoV- 
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NC99 strain was propagated in HRT-18G cells to make virus stocks 
for MN assay. ECoV MN assays were carried out by mixing 100 pi of 
serial 2-fold dilutions (1:10 to 1:640) of heat-inactivated sera (56°C for 
30 min) with equal volumes of 200 tissue culture infectious doses 
(TCID 50 ) of ECoV in 96-well microtitre plates (Cat:92096, TPP tissue 
culture plates, Switzerland). Plates were incubated for 1 hr at 37°C. 
The virus-serum mixture was then added in quadruplicate to HRT- 
18G cell monolayers in 96-well microtitre plates followed by 1 hr 
adsorption at 37°C. The plates were incubated at 37°C in 5% C0 2 in a 
humidified incubator for 5 days, when virus back titrations were 
inspected to confirm that the virus dose used was within the expected 
range. The antibody titre was the highest serum dilution that resulted 
in complete neutralization of cytopathic effect in >2 of the four wells 
replicate wells at each serum dilution. Antibody titre refers to the dilu¬ 
tion of the serum in the serum dilution-virus mixture. All antibody 
titres are denoted as the reciprocal of the titre. 

BCoV is ubiquitous in cattle worldwide and is genetically and 
antigenically related with HKU23, a coronavirus endemic in dromed¬ 
aries (Woo et al., 2014, 2016). Antibody to BCoV was tested using 
MN assay as described previously (Park et al., 2015; Perera et al., 
2013). Serology for BCoV and ECoV was carried out with culture 
medium without foetal calf serum. Known positive and negative con¬ 
trol sera were included in all serology assays. 


2.6 | Reference sera used in investigation of cross¬ 
reactivity in serological assays 

Cross-reactive responses in the MERS-CoV ppNT, BCoV MN and 
ECoV MN assays were assessed using a panel of previously charac¬ 
terized immune or hyperimmune reference antisera generated 
against alpha, beta and gamma coronaviruses. These sera were either 
obtained from BEI Resources (animal CoV reagents supplied to BEI 
by Dr Linda Saif http://www.beiresources. 0 rg/About/BEIResources.a 
spx) or kindly provided by Dr. Linda Saif or Dr Stanley Perlman. 
Equine sera known to be seropositive to ECoV were provided by Dr. 
Nicola Pusterla, Department of Veterinary Medicine and Epidemiol¬ 
ogy, School of Veterinary Medicine, University of California, Davis, 
CA, USA. The homologous antibody titres to the immunized virus 


were also obtained, where available, from the providers of these 
respective antisera. The serum panel for alpha coronavirus included 
sera generated against porcine respiratory coronavirus, feline infec¬ 
tious peritonitis virus, canine enteric coronavirus and porcine trans¬ 
missible gastroenteritis virus. Sera produced against mouse hepatitis 
virus strains JHM and A59, SARS coronavirus, BCoV, ECoV repre¬ 
sented the betacoronavirus group. Furthermore, an immune serum 
for infectious bronchitis virus was used as representative for the 
gammacoronavirus group. All these reference sera, as well as addi¬ 
tional ECoV immune sera and MERS-CoV-negative camel sera, were 
tested in ppNT or MN assays in parallel against MERS-CoV, BCoV, 
as well as ECoV, to assess possible cross-reactive antibody 
responses among these three coronaviruses. Some of these refer¬ 
ence sera, with the exception of horse and camel sera, had been 
previously tested against MERS-CoV and BCoV, and the results were 
reported by us (Hemida et al., 2014). 


3 


RESULTS AND DISCUSSION 


All nasal and rectal swabs were tested using the pan-coronavirus 
nested RT-PCR assay, ECoV and MERS-CoV upE gene RT-qPCR 
assays. All swabs were negative in the MERS-CoV RTqPCR assay. 
Five rectal swabs (5/315) were positive for coronavirus RNA by the 
pan-coronavirus assay (Tables 1 and 2). Sequencing of the 440 base 
pair nucleotide fragment from the RdRp gene confirmed four of 
these to be ECoV and the other to be dromedary camel HKU23-like 
virus designated as equine HKU23 CoV SA1354. Additional 
sequence was obtained from these specimens to generate a nucleo¬ 
tide sequence of 933 nucleotides (position 14,780 to 15,712 of 
ECoV strain NC99) which was used for analysis of nucleotide differ¬ 
ences in comparison with other coronaviruses (GenBank accession 
numbers KY458228-KY458232) (Table 3) and phylogenetic analysis 
which confirmed the identification of the coronaviruses to be ECoV 
and HKU23, respectively (Figure 1). This result led to the testing of 
306 nasal and 315 rectal swab specimens by specific RTqPCR assays 
for ECoV and RTqPCR assay for HKU23. Testing by the ECoV speci¬ 
fic RTqPCR led to the detection of five rectal swabs specimens, four 


TABLE 2 RT-qPCR-positive specimens 


Swab type 

Location and 
provider code 

Pan CoV nested 

RT-PCR result 
with sequencing 
of PCR product 

MERS-CoV 

UpE_RT-qPCR 

EqCoV_RT- 

qPCR.Ct 

HKU23 RT¬ 
qPCR Ct value 

HKU23 RT¬ 
qPCR Tm a 

Nasal 

Al-Quassem (10) 

Negative 

Negative 

Negative 

37.46 

82.26 

Rectal 

Dammam (3) 

HKU23 

Negative 

Negative 

34.67 

82.25 

Rectal 

Riyadh (1) 

Equine CoV 

Negative 

35.99 

Negative 


Rectal 

Riyadh (1) 

Equine CoV 

Negative 

29.38 

Negative 


Rectal 

Riyadh (1) 

Equine CoV 

Negative 

30.63 

Negative 


Rectal 

Riyadh (1) 

Negative 

Negative 

36.32 

Negative 


Rectal 

Riyadh (1) 

Equine CoV 

Negative 

34.87 

Negative 



a Tm of camel HKU23 is 82.138 degree C. Screening horse samples for HKU23 RTqPCR allowed Tm range 81.638 to 82.638 degree C. 










HEM I DA et al 


Jransboundary and Emerging Diseases 


Wiley 



Ov 


00 o 

O CN 


O 00 o 
odd 


I I 

o o o on 


(N rH H ^ d 

o o o o c\i 


CD 


> 

> 

> 

> 


o 

o 

o 

o 

CN 

O 

O 

O 

U 

CN 

CD 

CD 

CD 

CD 

U 

C 

C 

C 

.c 

Z 

3 

3 

3 

b 

s 

cy 

CT 

CT 

.CT 


LU 

LU 

LU 

LU 

o 

u 

O 

CN 

NJ- 

CN 


'■O 

ND 

"O 

00 

CD 

00 

00 

00 

00 

C 

T-H 

r-H 

r-H 

r-H 

o 

<C 

< 

< 

<C 

CT 

to 

to 

to 

to 

LU 


00 

CN 

r> 

"w 

£ 

N- 

U0 

00 


s 


ON 

O 


ON 


ON 

ON 


ON 

ON 


00 

ON 


00 

ON 


un 

NO 

ON 

00 

ON 

D 

X 

> 

o 

U 

ID 

E 

03 

U 










1 

LO 









T”i 








1 

Nj- 








NO 

o 








NO 

00 







1 

tH 

Nf- 







LO 

00 

ON 







T—1 

o 

00 






1 

Nf 

NT 

N]- 






hv 


NO 

q 






o 

ON 

K 

T—1 





1 

CN 

Nf 

00 






q 

o 

ON 

Dv 







CN 


ON 

ON 




1 

CN 



nT 

■sT 




t-H 

NO 

ni- 

''d' 

q 

NO 




CN 

CN 

o 

00 

NT 

00 



1 

Nf 

00 

nT 



"sT 



00 

CN 

CN 


00 

q 

00 



ON 

ON 

00 

t-H 

00 

LO 

00 


1 


Nt" 

00 

nt 


T—1 



°o 

LO 

00 

NO 

1 

ON 

q 

NO 


tH 

00 

ON 

CN 

tH 

Nf" 

LO 

00 

1 



Nt - 

00 

NT 

t-H 

T—1 


CN 

°o 

ON 

00 

LO 

tH 

o 

NT 

"sT 

o 


00 

CN 

CN 


NT 

LO 

00 




NT 

00 

"st" 

t-H 

t-H 

”sT 

o 

CO 

LO 

00 

NO 


CN 

q 

NO 

o 


00 

ON 

CN 


NT 

LO 

00 




NT 

00 

nT 



"sT 

CN 

NT 

00 

LO 

t-h1 

nT 

o 

cq 

NO 

t-H 


00 

CN 

CN 

o 

NT 

LO 

00 




NT 

00 

NT 

t-H 


"sT 

O 

un 

LO 

O 

00 

NO 

q 

q 

”sT 

ON 

d 

00 

00 

CN 

o 

NT 

LO 

00 




NT 

00 

nT 


t-H 

"sT 

O 

un 

LO 

O 

00 

NO 

q 

q 


ON 

t-H 

00 

00 

CN 

o 


LO 

00 




Nf 

00 

nT 


vH 

"sT 


hv 

NO 

CN 

LO 

00 

q 

t-H 

■sT 

ON 

T—1 

00 

00 

CN 

o 

NT 

NO 

00 




NT 

00 

Nt- 



NI¬ 


q 

NO 

ON 

LO 

00 

q 

q 

NO 

ON 

t-H 

00 

00 

CN 

o 

NT 

LO 

00 




nT 

00 

nT 

t-H 


nT 



rv 







00 


o 







CN 

rH 

\ 

1- 

d 




X 



DC 

Z 

NO 




tin 



LU 

LO 




D 



00 

CN 

Z) 

X 

> 

O 

u 

CD 

00 

2 

o 

to 

QC 

LU 

nj- 

z> 

LO 

D 

L_ 

*> 

<S) 

+-» 

D 

X 

to 

QC 

< 

to 


> 

> 

> 



03 

> 

> 

> 

o 

o 

O 

X 

X 

L_L 

CD 

J= 

o 

o 

o 

U 

U 

u 

U 

> 

> 

u 

u 


CD 

c 

c 

o 

o 

CD 

c 

c 

"a; 

C 

CD 

03 

(J 

(J 

in 

(33 

(33 

E 

03 

> 

o 

E 

Z3 

E 

13 

+-> 

(33 

(33 

D 

o 

E 

E 

Z3 

U 

CD 

X 

X 

CD 

CD 


X 

X 








HEM I DA et al. 


Wiley 


■>fransboundary and Emerging Diseases 


1 \ 


-HCoVHKUl (NC 006577) 

Equine CoV SA1360 
Equine CoV SA1362 
Equine CoV SA1364 
Equine CoV SA1382 
l~ Equine CoV strain NC99 (EF446615) 

_ Equine CoV RNA strain Tokachi09 (LC061272) 

I Equine CoV RNA strain Obihirol 2-2 (LC061274) 

Equine CoV RNA strain Obihirol 2-1 (LC061273) 

— Porcine hemagglutinating encephalomyelitis virus strain VW572 (DQ011855) 

HCoV 0C43 (AY391777) 

HCoV 0C43 serotype 0C43-Paris (AY585229) 

HCoV 0C43 strain ATCC VR-759 (AY585228) 

HCoV 0C43/human/USA/851-15/1985 8-Jan-85 USA Nashville TN (KF530060) 

HCoV 0C43/human/USA/873-16/1987 12-Mar-87 USA Nashville TN (KF530077) 

HCoV 0C43/human/USA/873-19/1987 17-Mar-87 USA Nashville TN (KF530083) 

HCoV 0C43/human/USA/872-5/1987 10-Feb-87 USA Nashville TN (KF530086) 

Canine respiratory CoV strain K37 2008 lung South Korea (JX860640) 

Bovine CoV genomic RNA strain Kakegawa (AB354579) 

Bovine CoV strain Quebec (AF220295) 

Bovine CoV strain Mebus (U00735) 

Human enteric CoV 4408 1988 Germany (FJ415324) 

Human enteric CoV strain 4408 17-Apr-09 (FJ938067) 
r BetaCoV 1 isolate Buffalo CoV B1-24F 1-Oct-14 Bangladesh (KU558922) 

BetaCoV 1 isolate Buffalo CoV B1-28F 1-0ct-14 Bangladesh (KU558923) 

Dromedary camel CoV HKU23 strain HKU23-265F 3-Apr-13 UAE Dubai (KF906249) 
Equine HKU23 CoVSA1354 

Dromedary camel CoV HKU23 strain HKU23-368F 21-May-13 UAE Dubai (KF906251) 
Bovine CoV E-DB2-TC 30-Nov-96 USA Ohio (FJ938063) 

Bovine CoV DB2 1983 (DQ811784) 

~ Bovine CoV strain BCV-AKS-01 30-May-15 China Xinjiang (KU886219) 

White-tailed deer CoV US/0H-WD470/1994 1-Jan-94 USA Ohio (FJ425187) 

Bovine respiratory CoV bovine/US/0H-440-TC/1996 1-Jan-96 USA Ohio (FJ938066) 
Sambar deer CoV US/0H-WD388/1994 1-Jan-94 USA Ohio (FJ425189) 

Sambar deer CoV US/0H-WD388-TC/1994 1-Jan-94 USA Ohio (FJ425188) 

Waterbuck CoV US/0H-WD358-GnC/1994 1-Jan-94 USA Ohio (FJ425185) 

Waterbuck CoV US/0H-WD358-TC/1994 1-Jan-94 USA Ohio (FJ425184) 

Waterbuck CoV US/0H-WD358/1994 1-Jan-94 USA Ohio (FJ425186) 


0.09 


Equine coronavirus 


Human coronavirus OC43 


Bovine coronavirus/HKU23 cluster 


FIGURE 1 Phylogenetic analysis of equine coronavirus, human coronavirus OC43 and bovine coronavirus/dromedary camel coronavirus 
HKU23 cluster. The maximum likelihood phylogenetic tree rooted to HKU1 as an outgroup was generated using virus RdRp gene alignment of 
934 nt by PhyML with approximate likelihood ratio test (aLRT) SH-like branch support Sequences generated from this study are highlighted in 
red. GenBank accession numbers of are indicated in brackets. 


TABLE 4 MERS-CoV serology results 

, Reciprocal antibody titre 

Location and = _ _ 3 _ 


Horse ID 

provider code 

MERS-ppNT 

MERS MN 

MERS PRNT90 

MERS PRNT50 

BCoV MN 

ECoV MN 

5-36 

Al-Hasa (5) 

40 

<20 

<20 

20 

<20 

80 

5-41 

Al-Hasa (5) 

40 

<20 

<20 

20 

<20 

20 

9-2 

Al-Qassem (9) 

160 

20 

40 

80 

<20 

<20 

13-12 

Oman (13) 

40 

<20 

<20 

20 

<20 

<20 

2-12 

Riyadh (2) 

40 

<20 

<20 

20 

20 

40 


ppNT, pseudotype neutralization test; MN, microneutralization test; PRNT, plaque reduction neutralization test. 


of which had been detected by the pan-coronavirus RT-PCR assay 
as well as an additional positive rectal swab (Ct value 36.3) from a 
horse in the Riyadh area (Table 2). Using the HKU23-specific 
RTqPCR, two positives were detected including SA1354 identified as 
HKU23 by sequencing the pan-coronavirus RT-PCR assay amplicon 
and an additional positive result with a high Ct value of 37.6 from a 
nasal swab sample. None of the other nasal swabs had evidence of 
coronavirus RNA. Compared to the pan-coronavirus assay, the speci¬ 
fic ECoV- and HKU23-specific assays detected one additional speci¬ 
men that was weakly positive for each of these viruses (Table 2). 

The five pan-coronavirus PCR CoV-positive specimens had Ct 
ranging from 29.4 to 35.9 in either the ECoV or HKU23 specific 
assays. The pan-coronavirus RT-PCR was not sensitive at detecting 
specimens with Ct higher than 36 in the ECoV or HKU23 assays. 


However, the results do demonstrate the importance of using the 
pan-coronavirus assay when investigating coronavirus activity in ani¬ 
mal species in new geographic areas. Our initial focus was to explore 
for evidence of MERS-CoV infections in horses in Saudi Arabia and 
Oman. None of the rectal swabs were positive for MERS-CoV RNA, 
but the application of the pan-coronavirus RT-PCR to these speci¬ 
mens revealed the presence of two other coronaviruses infecting 
horses in this region, that is ECoV and HKU23. 

All seven horses from whom ECoV or HKU23 were detected 
were apparently healthy. The specimens positive for ECoV RNA all 
came from different horses at a horse breeding facility in Riyadh. 
The two horses positive for HKU23 RNA were from Dammam and 
Qateef, Saudi Arabia, where horses were used for farming, and in 
Dammam, also for breeding (Tables 1 and 2). In both cases, these 
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TABLE 5 Cross-neutralization titres (denoted as reciprocal titres) for Middle East respiratory coronavirus (MERS-CoV), bovine coronavirus 
(BCoV) and equine coronavirus (ECoV) in hyperimmune or naturally infected sera known to be positive for different coronaviruses 


Genus 

Sera panel—Catalogue numbers provided 
for sera obtained from BEI Resource 

Homologous antibody titre by 
ELISA unless specified 

MERS- 

ppNT 

MERS 

MN 

BCoV 

MN 

ECoV 

MN 

Alpha coronavirus 

NR2727 - Guina pig antiserum to canine 

coronavirus 

4,094 b 

<20 

<20 

<20 

<20 


NR460 - pig antiserum to porcine respiratory 
coronavirus 

l,200 a 

<20 

<20 

<20 

<20 


NR2518 - Guina pig antisera feline infectious 
peritonitis virus 

2,000 a 

<20 

<20 

<20 

<20 


NR458 - pig antisera for porcine transmissible 
gastroenteritis virus 

l,400 a 

<20 

<20 

<20 

<20 

Beta coronavirus 

ECoV-positive sera 459 (horse, USA) 

qPCR negative; ELISA OD 1.459 d 

<20 

<20 

160 

>1,280 


ECoV-positive sera 463 (horse, USA) 

qPCR positive; ELISA OD 1.525 d 

<20 

<20 

1,280 

>1,280 


ECoV-positive sera 472 (horse, USA) 

qPCR negative; ELISA OD 1.458 d 

20 

<20 

640 

>1,280 


ECoV-positive sera 486 (horse, USA) 

qPCR negative; ELISA OD 1.939 d 

20 

<20 

1,280 

>1,280 


ECoV-negative sera 2 (horse, Hong Kong) 

Not relevant 

<20 

<20 

<20 

<20 


BCoV-positive camel serum #740293 e 

Not relevant 

<20 

ND 

640 

80 


BCoV-positive camel serum #467468 e 

Not relevant 

<20 

ND 

80 

<20 


C99-10 - BCoV antisera from guinea pig 

20,480 b 

<20 

<20 

320 

<20 


Bcov antisera from germ free bovine calf 

580 b - neutralisation titre 

<20 

<20 

1,280 

80 


NR456 - BCoV antisera from gnotobiotic calf 

10,000 a 

<20 

<20 

80 

<20 


NRC772 - Rabbit antisera SARS S protein high titre 

640 

<20 

<20 

<20 

<20 


Mouse hepatitis virus(JHM strain) hyperimmunized 
mouse dam 1 

l,778 c - neutralisation titre 

<20 

<20 

<20 

<20 


Mouse hepatitis virus(JHM strain) hyperimmunized 
mouse dam 2 

363 c - neutralisation titre 

<20 

<20 

<20 

<20 


Mouse hepatitis virus(A59 strain)-infected mice 

l,000 c - neutralization titre 

<20 

<20 

<20 

<20 


NRC774 - antisera for SARS coronavirus sero titre 

<10 

<20 

<20 

<20 

<20 


NRC769 - Rabbit antisera to SARS S protein sero 
titre 

<10 

<20 

<20 

<20 

<20 

Gamma coronavirus 

NR2515-Guina pig antiserum to infectious 
bronchitis virus 

50,000 a 

<20 

<20 

<20 

<20 


ppNT, Pseudoparticle neutralization test; MN, micro neutralization test; PRNT, plaque reduction neutralization test; BCoV, bovine coronavirus; ECoV, 
equine coronavirus; OD, optical density. 

Equine antisera against ECoV were kindly provided by Professor Nicola Pusterla, Department of Veterinary Medicine and Epidemiology, School of 
Veterinary Medicine, University of California, Davis, CA, USA. 

Homologous titres or OD results obtained from BEI resources, Linda Saif, Stanly Perlman and Udeni Balasuriya represented as a ’ b> c> d respectively. 
e Field collected dromedary sera from Kazakhstan, a region known to be free of MERS-CoV infection. 


horses frequently came into contact with camels and other domestic 
livestock. HKU23 coronavirus is known to be endemic in camels 
(Woo et al., 2014). 

We investigated for serologic evidence of MERS-CoV infections 
in the 243 horse sera using ppNT screening followed by endpoint 
titrations. Any sera positive at a screening reciprocal antibody titre 
of 20 in ppNT were confirmed using MN assay and PRNT assay 
using live MERS-CoV strain EMC. Possible serologic cross-reactions 
with ECoV and BCoV were evaluated in MN assays which are the 
most specific serological assays available. Of the 243 sera tested, 
only five (2%) had detectable MERS-CoV antibody in the MERS- 
ppNT assay at a screening titre of 20; four of these were from Saudi 
Arabia, and one was from Oman. The endpoint MERS-ppNT titres of 
these five sera ranged from 40 to 160 but only serum 9-2 which 


had the highest MERS-ppNT titre was positive in the MERS-CoV 
MN test (titre 20) (Table 4). In the confirmatory PRNT which is con¬ 
sidered the “gold standard” serological test for MERS-CoV, horse 
serum 9-2 reduced plaque numbers by 90% (PRNT 90 ) up to a titre 
of 40. This serum was negative in ECoV and BCoV MN assays 
(Table 4). This serum was from a racing horse stable in Al-Qassem 
and may represent a rare example of transmission of MERS-CoV to 
horses or cross-reaction with another yet undocumented coronavirus 
infecting horses. This horse was not known to frequently come into 
contact with camels. Three of the other sera (serum ID nos 5.35; 
5.41; 2.12) positive in the MERS-CoV ppNT assay had detectable 
ECoV MN antibody and one of these also had detectable BCoV MN 
antibody (Table 4). Thus, these may well represent evidence of 
cross-reactive responses. 
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TABLE 6 Seropositivity to BCoV and ECoV in horses from different geographic locations 


Location 

Provider code and 

management system 

Sample ID 

Reciprocal 

BCoV MN titre 

Reciprocal 

ECoV MN Titre 

PCR results 

Riyadh 

1 (Breeding) 

l-l a 

<20 

20 

ECoV PCR positive 



l-2 a 

<20 

20 

ECoV PCR positive 



vl-3 a 

20 

160 

ECoV PCR positive 



1-4 

<20 

80 




1-5 

<20 

160 




1-6 

<20 

40 




l-8 a 

<20 

640 

ECoV PCR positive 



l-12 a 

80 

320 

ECoV PCR positive 


2 (Farming, breeding) 

2-13 

40 

160 




2-14 

<20 

80 


Dammam 

3 (Farming, breeding) 

3-11 

<20 

<20 




3-12 

<20 

40 




3-41 a 

<20 

<20 

HKU23 PCR positive 

Al-Hasa 

5 (Breeding) 

5-56 

<20 

40 




5-57 

160 

320 




5-58 

40 

80 




5-59 

80 

160 



8 (Farming, breeding) 

8-2 

20 

80 




8-3 

1:80 

40 




8-4 

<20 

640 




8-5 

<20 

<20 




8-6 

<20 

20 



7 (Farming) 

7-1 

<20 

20 




7-2 

20 

<20 




7-3 

<20 

<20 




7-4 

<20 

20 




7-1 

<20 

<20 




7-2 

<20 

<20 


Al-Qassem 

9 (Racing) 

9-A3 

20 

<20 




9-A4 

<20 

40 




9-A5 

<20 

<20 




9-B1 

<20 

<20 




9-C1 

20 

20 




9-C2 

80 

640 




9-El 

<20 

<20 




9-E2 

<20 

<20 


Qateef 

10 (Farming) 

10-13 

<20 

80 




10-14 

<20 

20 




10-15 

<20 

80 




10-16 

20 

160 



11 (Farming) 

11-8 

80 

40 




11-9 

<20 

320 




11-11 

<20 

<20 




11-12 

<20 

<20 




11-13 

<20 

80 



(Continues) 
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TABLE 6 

(Continued) 





Location 

Provider code and 

management system 

Sample ID 

Reciprocal 

BCoV MN titre 

Reciprocal 

ECoV MN Titre 

PCR results 

Aljobail 

12 (Farming) 

12-1 

40 

40 


Oman 

13 (Racing, breeding) 

13-1 

<20 

20 




13-2 

<20 

80 




13-3 

<20 

40 




13-4 

20 

160 




13-5 

20 

40 




13-6 

20 

160 




13-10 

20 

40 




13-11 

<20 

80 



’Horses found to be coronavirus RNA positive in swabs by RT-PCR or RT-qPCR. 


The extent of the serological cross-reactivity of MERS-CoV, 
ECoV and BCoV was systematically investigated. We did not have 
an isolate of HKU23 virus available for serological testing and BCoV 
was used as a virus that is genetically closely related to HKU23 with 
pairwise amino acid similarity of 94.1% in the spike protein gene, 
the determinant of antigenic cross-reaction in neutralization tests 
(Woo et al., 2014). Horse sera obtained from the USA reported to 
be strongly reactive by ELISA (optical density ranging from 1.458 to 
1.939) for ECoV showed high titre (>1,280) for ECoV in MN assay 
as expected, but also had high antibody titres (ranging from 160 to 
1,280) in BCoV MN tests (Table 5). As these horse sera were from 
horses in the “field,” they may well have been exposed to multiple 
coronaviruses. Thus, these results may reflect cross-reactivity or co- 
infection of US horses with ECoV- and BCoV-like viruses. Two of 
these sera showed low titres (20) in the MERS-CoV ppNT assay but 
were negative in the more stringent MERS-CoV MN assay. As 
MERS-CoV is not circulating in the USA, these results suggest cross¬ 
reactivity between MERS-CoV and ECoV or related CoVs may occur, 
albeit at low titre. Similarly, a BCoV immune serum from an experi¬ 
mentally infected gnotobiotic calf showed detectable, but 16-fold 
reduced antibody titre with ECoV but no cross-reaction with MERS- 
CoV. Spike proteins of equine coronavirus and bovine coronavirus 
have an overall amino acid similarity of 80.8% and 73.1% similarity 
of the SI region which contains the receptor-binding domain and 
neutralization epitopes. This probably explains the large difference in 
neutralization titres observed in Table 5. A dromedary serum from 
Kazakhstan, a region known to be free of MERS-CoV activity 
(Miguel et al., 2016), had detectable antibody to BCoV (probably 
reflecting HKU23 or BCoV infection) and also to ECoV. We cannot 
conclude whether this reflects cross-reactivity or infection with 
BCoV, HKU23 and/or ECoV. None of the beta coronavirus immune 
sera to SARS-CoV or mouse hepatitis virus, nor immune sera to 
alpha or gamma coronaviruses gave cross-reactions in BCoV or 
ECoV MN assays (Table 5). As had been previously reported, none 
of these dromedary sera cross-reacted with MERS-CoV (Hemida 
et al., 2014). 

Fifty-four horse sera from Saudi Arabia and Oman were selected 
to represent different collection locations, excluding those five sera 


previously reported to be MERS-CoV ppNT positive (Table 6). These 
sera were seronegative for MERS-CoV (as expected); 40 (74%) of 
them had detectable MN antibody to ECoV (titres ranging from <20 
to 640; median 40), and 18 (33.3%) had MN antibodies to BCoV 
(titres ranging from <20 to 160; median <20). The reciprocal geomet¬ 
ric mean antibody titres (assigning sera with titres <20 the value of 
10 for computational purposes) for ECoV and BCoV were 44 and 
25.5, respectively. There were twenty-four sera with undetectable 
antibody to BCoV with detectable ECoV antibody titres ranging from 
20 to 640. The titres to ECoV were higher than or equal to titres to 
BCoV, with two exceptions where the titre to BCoV (20) was mar¬ 
ginally higher than that with ECoV (Table 6). These results are com¬ 
patible with the RT-PCR detection of ECoV in rectal swabs in horses 
and probably suggest that ECoV or a closely related virus is com¬ 
monly infecting horses in many regions of the Arabian Peninsula. 
Two of the five horses in whom ECoV RNA was detected by 
RTqPCR had low ECoV MN antibody titres (20) while the other 
three had titres ranging from 160 to 640 (Table 6). This may reflect 
virus shedding persists into convalescence or that re-infection of 
previously seropositive horses can occur. 

The RTqPCR detection of two HKU23-like viruses indicates that 
HKU23 which is common in dromedaries in the region is also infect¬ 
ing horses. In the absence of an HKU23 isolate for use in neutraliza¬ 
tion tests and given the potential serological cross-reactivity 
between these three viruses, it is difficult to conclude how com¬ 
monly equines are infected with HKU23. The antibody titres to 
BCoV, which is genetically closely related to HKU23 in the spike 
protein gene, may suggest that BCoV- or HKU23-like virus infection 
of horses does occur in Saudi Arabia. 

The serological evidence of one of the 243 sera with detectable 
MERS-CoV antibody titres in the PRNT 90 test raises the question 
whether MERS-CoV may occasionally cross species to horses. This 
serum did not have any detectable antibody to ECoV or BCoV. 
However, it is not possible to exclude the possibility that this obser¬ 
vation reflects cross-reaction with ECoV, BCoV or other yet 
unknown coronaviruses. Confirmation of natural infection of horses 
with MERS-CoV can only come from detection of MERS-CoV RNA 
in horse specimens. Recent experimental infection of horses with 
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MERS-CoV did not result in virus replication or seroconversion 
(Adney et al., 2016). 

In conclusion, the serological data suggest that lineage A 
betacoronaviruses are commonly infecting horses in Saudi Arabia 
and Oman but antibody cross-reactions between these viruses do 
not permit us to use serological data alone to identify which 

coronaviruses are causing these infections. RT-PCR detection of 
ECoV and HKU23 in equine swabs confirms the circulation of 

these two viruses in horses in Saudi Arabia. Cocirculation of 

related viruses may provide potential for recombination, a poten¬ 
tial means of generating genetic diversity and facilitating host 

jumps in coronaviruses. This is the first report of HKU23 being 
detected in horses and the first detection of ECoV in Asia outside 
of Japan. 
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